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Abstract 



We study jet final states in diffractive deep-inelastic scattering. We present 
the QCD factorization formula in terms of diffractive parton distributions 
and discuss its implementation in NLO Monte Carlo generators. We compute 
NLO predictions for the diffractive jet cross section. We use this calculational 
framework to discuss theoretical models of the long-distance behavior. 
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1. Introduction 



Among the most intriguing results that have come from the HERA lepton- 
hadron colhder are observations of diffraction hard scattering [f]. In these 
processes, while the short-distance behavior is the one well known from or- 
dinary deep inelastic scattering (DIS), the long-distance matrix elements are 
much more complex than the usual parton distributions [2, 3] and fully in- 
corporate the dynamics of hadron diffraction. 

These diffractivc matrix elements, or diffractive parton distributions, have 
so far been considered mostly in the context of structure functions. One of 
their key properties, however, concerns the size of the gluon distribution, 
which, in sharp contrast with the case of ordinary DIS, dominates the quark 
distributions by one order of magnitude even at low mass scales and large 
momentum fractions. This enhances the importance of processes, such as the 
production of jets, that unlike structure functions couple directly to gluons. 

The purpose of this paper is to start to analyze systematically the jet 
structure of diffractive final states in terms of diffractive parton distributions 
and ncxt-to-leading-order (NLO) short-distance cross sections. A generic jet 
observable a will be represented schematically as 



where f^^"^ is the diffractive parton distribution, H is the hard scattering 
function, calculable as a power series expansion in a^, and W is the mea- 
surement function, specifying the definition of the observable in terms of the 
final-state kinematic variables. QCD analyses of diffractive jet production 
have so far been limited to the leading logarithms [4, 5, 6, 7]. In this paper, we 
present the QCD factorization formula that allows nonleading corrections to 
be systematically included and we carry out the calculation at the NLO. We 
find that contributions beyond the leading logs are important in the HERA 
kinematic region. 

Wc argue that, once pcrturbativc corrections arc consistently taken into 
account, jet production can be used to test theoretical models of the long- 
distance behavior of diffraction scattering. Motivated by a previous com- 
parison [8] of theory with diffractive F2 data [9, 10], we consider a simple 
physical picture [3, 8] in which the diffractive gluon distribution results pre- 
dominantly from color-octet dipoles penetrating the target at small trans- 
verse separations of order !/«;, where 1/k represents a color transparency 
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length. This picture gives testable predictions, because the smallness oil/ k 
justifies a perturbative calculation of the j3 dependence of the distributions. 
We compute predictions for jet cross sections in this scenario. 

The paper is organized as follows. We start in Sec. 2 by presenting the 
factorization formula for diffractive jet production and its implementation 
in NLO Monte Carlo generators. This discussion provides a "^-channel" 
picture of the process, based on the separation of high and low transverse 
momenta at scale /i, and renormalization-group evolution in yU. In Sec. 3 we 
examine the diffractive parton distributions. To gain insight into the form 
of these distributions, it is useful to exploit the intuition that comes from an 
"s-channel" picture of the process, based on the target rest frame and the 
light-cone evolution of the parton system created by the operator that defines 
the distribution. In Sec. 4 we use the results of Sees. 2 and 3 to compute the 
diffractive jet cross section, and present numerical results at the NLO. Due to 
our lack of familiarity with the kinematic cuts employed in the experiments, 
we leave to the experimental groups the comparison of these results with the 
data. We summarize and give concluding remarks in Sec. 5. 



2. NLO event generators and hard diffraction 

In this section we discuss the diffractive factorization formula. Using 
properties of this formula, we adapt standard NLO Monte Carlo generators 
for DIS in order to perform NLO calculations of diffractive hard processes. 

Consider the diffractive leptoproduction of jets, e + A^e' + A' + J + X, 
by virtual photon exchange. We parameterize the momenta of the incoming 
hadron and virtual photon in a coUinear reference frame as 

where is the photon virtuality, x is the Bjorken variable, and we have used 
lightcone components = ±p^)/v^- The Breit frame, in which we will 
define the jets, is identified by p\ = Q/{x\/2). We set y = Q'^/{xS), with S 
the lepton-hadron center-of-mass energy squared. Hadron A' is characterized 
by the fractional loss of longitudinal momentum Xp — 1 —Va'IVa 
invariant momentum transfer t — (p^ — p^/)^. 
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To be definite, we suppose fixing Q'^ and y in the Icpton subprocess, Xp 
in tiic difTractive subprocess, and the transverse energy of the jet, whose 
precise definition is specified below (see Sec. 4). We integrate over t up to a 
maximum value tmax, which we take to be of the order of a GeV. We consider 
the fourfold-differential jet cross section da /[dETdQ'^dydxp]. According to 
the hard-scattering factorization theorem [2, 3, 11, 12], this cross section is 
given, up to corrections suppressed by powers of the hard-scattering scale, 

by 



--— -— (ET,(3^5,2/,a;f) = J2 f 
dErdQ^dydx^ a^g,gJp P' 



where (3 is given by 

dfa^ydxp is the diffractive parton distribution integrated over t, and d&a/ 
[dETdQ'^dy] is the partonic jet cross section, incorporating the hard-scattering 
function H of Eq. (1) (through perturbatively calculable matrix elements) 
and the measurement function W (through the algorithm that defines the 
jet). The mass scales /ir and /i are the renormalization and factorization 
scales. 

Eq. (3) expresses the separation of short- distance and long-distance con- 
tributions to the diffractive jet cross section. It contains two crucial dif- 
ferences compared to the analogous factorization formula for inclusive jet 
production: a) the distributions f^^\ which embody all of the dynamical 
effects of hadron diffraction, and b) the diffractive kinematics, which enters, 
through the momentum loss Xp, not only in f^^^ but also in (3, Eq. (4). The 
definition of f^^^ can be found in [2, 3] and is briefly recalled in Sec. 3. 

The important point for practical applications is that the dependence 
on the hard scale, Q or Et, factors out of f^^^ in Eq. (3): it is entirely 
contained in 0"^, as in the inclusive case. Since the left hand side of Eq. (3) 
is independent of the factorization scale /x, this also implies that, as in the 
inclusive case, the evolution of /^^^ with /j, is given by renormalization- group 
equations. Because we do not integrate over large t (of the order of the hard 
scale or higher), the form of these equations is simply DGLAP, without the 
inhomogeneous term from extra large-i subtractions [2] . 
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We can thus calculate the differential cross section (3) to the next-to- 
leading order using the NLO results that are available for the partonic jet 
cross sections (Xa, and NLO evolution of the f^^\ The cross sections aa 
start at order a^. The order-a^ contributions were computed and encoded in 
Monte Carlo programs by three groups [13, 14, 15]. Two of these programs, 
DISASTER++ [13] and DISENT [14], have been found [16] to be in good 
numerical agreement. More recently two other independent programs [17, 18] 
have appeared. The authors of the program NLOJET++ [18] have also per- 
formed a comparison of different codes and found agreement of their results 
with those of [13, 14]. Any one of these Monte Carlo programs can be used 
to evaluate the short distance part of Eq. (3) to NLO ^. 

The numerical results that we present in this paper are obtained using 
DISENT. This program generates DIS events according to a structure rem- 
iniscent of parton shower algorithms [21]: first, it generates a parton-model 
event, with one parton in the final state; then, starting from this and using 
the dipole subtraction method [14], it generates events with two partons in 
the final state; then, in a similar way it generates events with three final par- 
tons. Contributions to the jet cross section come from events with at least 
two partons in the final state. For applications to diffractive production, we 
generate two-parton and three-parton events so that and (3 of the starting 
event are preserved. Note that if we use variables Q^, y and P for generating 
the DIS phase space (any two of which can be fixed), the dependence on the 
diffractive momentum loss Xf, decouples from the convolutions implied by 
factorization and evolution. We will use this method for the calculations of 
Sec. 4. 

In the next section we turn to the long-distance part of the process. 

3. Diffractive parton distributions 

The diffractive parton distributions describe the nonperturbative dynam- 
ics of the hadronic state, and are defined in [2, 3] in terms of matrix elements 

^ Note that a numerical discrepancy between DISASTER++ and DISENT has been 
observed [19] in the case of certain event shape distributions, e.g. the jet broadening, and 
that DISASTER++ has been identified [20] as the program giving the correct result. As 
far as we can judge, this discrepancy does not affect the cross sections and the range in 
transverse energies considered in this paper. 
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of nonlocal quark and gluon field operators. For gluons one has 

-4^((3, x^, = j—f- V / dt dy-e'^^^^iy 

dxj, '^^ {Anypxr.p^ Jo J ^ 

x(A|(5„(0)+^|A',X)(A',X|G,(0,|/-,0)+^|A) , (5) 

with 

(5„(y)+^' = E{y) Q-^v{^-^9 f dx' AX{y-^ .x' ,y)t^ , 

(6) 

where G is the gluon field strength, A is the vector potential, t are the 
generators of the adjoint representation of SU{3), and is a color-octet 
eikonal line operator, with V denoting the path ordering of the exponential. 
An analogous definition applies for quarks. Unlike the ordinary (inclusive) 
parton distributions, the distributions /^^^ represent interactions that occur 
both long before and long after the hard scatter. 

The factorization of the previous section implies that the /*^^^'s can be 
measured from diffractive DIS data and used in Eq. (3) to predict jet cross 
sections. This provides an approach to diffractive jet production that is fully 
consistent, although agnostic about the form of the matrix elements (5), 
in the spirit of the parton model. Calculations based on this approach are 
carried out to leading order in [5, 6] and to NLO in this paper. From this 
point of view, experimental tests of Eq. (3) to NLO are an important goal of 
diffractive jet studies. 

A complementary approach consists of modeling the dynamics that de- 
termines the form of the matrix elements (5), and using Eq. (3) to test these 
models. It is this point of view that we now wish to take. To this end, 
we build on the observation [3, 22] that the main qualitative aspects of the 
data [9, 10] for the diffractive F2 structure function — steep rise with de- 
creasing Xp, fiat spectrum in (3, positive slope in up to /3 ~ 1/2 — are 
all consistent with the hypothesis that the diffractive gluon distribution be 
dominated by transverse momenta of the order of a (semi)hard scale Msh, 

MsH ~ 0(1 GeV) . (7) 

This scale has nothing to do with the hard-scattering scale, Q or Et, which 
has been factored out and does not appear in f^^\ Rather, it represents an 
intermediate scale between the scale of hard physics and Aqcdj whose origin 
is nonperturbative and associated with the hadron's soft color field. 
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The physical meaning of the scale Msh was explained in [8] in the refer- 
ence frame in which the target is at rest. In this frame, the scale M^n is seen 
to arise as a color-transparency scale, i.e., as the inverse of the maximum 
size 1/k for which the color-octet dipole system created by the operator (5) 
can go through the hadron without breaking it up. The analysis [8] finds 
that the transparency lengths are smaller for color-octet dipoles than for 
color-triplet dipoles. The hypothesis of semihard dominance may therefore 
be better verified for observables (such as jets) that couple directly to the 
gluon distribution than for structure functions. 

Note that the rest frame approach gives an "s-channel view" of the pro- 
cess, characteristic of color-dipole models. Note however an essential differ- 
ence between the treatment above and standard color-dipole models (see, 
e.g., Ref. [4] for a color-dipole study of diffractive jet production): in the 
present paper the s-channel picture is applied to the parton system created 
by the operator in f^^^ rather than to the system into which the virtual 
photon dissociates. That is, it is apphed after factorizing the hard-scattering 
subgraph. This difference is relevant, because the factorization provides a 
framework in which next-to-leading radiative corrections and evolution can 
be readily implemented. In contrast, such effects are difficult to include in 
color-dipole models, and evolution is typically taken into account at most to 
the leading- log level: see, e.g., [23] for a recent study. 

The semihard dominance picture leads to testable predictions, because the 
f3 dependence of the diffractive parton distributions then becomes calculable 
by a perturbation expansion. This calculation was first done in [22], using 
hght-cone perturbation theory. We observe that the semihard-scale scenario 
is common to several different models, including the large-nucleus model [24] 
and the saturation model [25]. In the following we use the distributions of 
[8], which incorporate the calculation [22]. At a starting mass scale ji — jiQ 
(with fiQ of order Msh), these distributions have the form 

1 1 

(8) 

where the ip^s are the perturbatively-computed functions, while Ag, Ag, a 
have to be determined from the data. The overall normalizations are pro- 
portional to the square of the color-transparency scales, Ag oc Hgr\ and 
Aq OC HqT^ [8], with ta the hadron radius. 

In Eq. (8) we have assumed, as is commonly done, a simple factorizing 
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form xlr^"' for the Xp dependence. This is an ansatz, not a theory result. It 
is in fact possible that the gluon and the quark distributions have distinct 
Xp behaviors at low fj, [26]. For illustrative purposes in the calculations that 
follow we limit ourselves to the simple factorizing form. 

Observe that the larger the larger the parton distribution: large f^^^ is 
related to the small size of the color-octet dipoles interacting diffractively [8] . 
The calculation of the functions and shows that the ratio of the gluon 
and quark distributions is proportional to a large color factor, C\{N^ — 
l)/{CpNc) = 27/2 [22]. In contrast with ordinary DIS, these functions do 
not have a fast fall-off as /3 — > 1, resulting in the gluon distribution being 
large even at large /3 [22]. 

As we have seen in Sec. 2, the evolution of the distributions (8) with the 
scale is given, up to power-like corrections, by the DGLAP equations, 
and is to be computed to the NLO, consistently with the accuracy of the 
partonic cross section. We do this using the moment-space evolution program 
of [27]. Numerical tables for the resulting distributions in the MS scheme are 
available from http:/ /zebu. uoregon.edu/~parton/diffpartons. 

The use of the DGLAP evolution equations is consistent with the leading 
power accuracy of Eq. (3). However, as emphasized earlier, the diffractive 
gluon distribution /^^^ is very large. Then corrections to factorization and 
evolution of relative order f^^^ /{taQY, corresponding to multi-parton cor- 
relations, although power-suppressed should likely be significant. Diffractive 
final states are therefore an especially important case in which to look for 
signals of behaviors beyond the leading power. 



4. Jet cross sections 

In this section we calculate the cross section for diffractive leptoproduc- 
tion of jets by evaluating Eq. (3) at the NLO. For the short- distance part 
of the process, we use the Monte Carlo program DISENT as discussed in 
Sec. 2. For the long-distance part, we use the diffractive parton distributions 
described in Sec. 3. 

We must specify the prescription for converting the final-state partons 
into jets. Beyond the leading logarithms, this becomes an essential ingredi- 
ent of Eq. (3): it is part of the function W in Eq. (1). Standard, infrared-safe 
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definitions are available both for cone jet algorithms and for /c^-clustering 
jet algorithms: see, e.g., [28] for a recent, comprehensive discussion. In the 
case of inclusive jet measurements [29] in DIS, the HI Collaboration have 
performed a comparison of different k± algorithms and found the algorithm 
of [30, 31] to have the smallest hadronization corrections. Since we do not 
expect any dramatic difference from the diffr active case in this respect, we fol- 
low this observation and present results for jets defined using this algorithm. 
As in [29], we adapt the algorithm, originally defined for hadron-hadron col- 
lisions, to the lepton-hadron case, working in the Breit frame. We use the 
Snowmass recombination prescription [32] to define the transverse energy Eti 
pseudorapidity r) and azimuth of the jet in terms of the momenta of its con- 
stituent partons. In what follows these kinematic variables are understood 
to be in the Breit frame. 
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Figure 1: The diffractive one-jet cross section at NLO as a function of the jet 
transverse energy Et for different bins and Xp — 0.05 (v^ = 300 GeV). 
The energy fraction y is integrated over the interval 0.1 < y < 0.7. Jets are 
defined by /c^-clustering, as specified in the text. 
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Figure 2: Size of the NLO correction. The values of the kinematic variables 
and range of integration are as in Fig. 1 (upper curve). 

In the calculations that we present below we integrate the momentum 
transfer t from t = to a maximum value t„iax = 1 GcV. We set the 
factorization and renormahzation scales to = A*|j = Q^- We set the jet 
resolution parameter R [30] to the value R = 1. 

In Fig. 1 we integrate the cross section (3) over y, with 0.1 < y < 0.7, 
and plot the triple-differential cross section da /[dETdQ'^dxp] versus Et at a 
fixed value of Xp for three different bins in Q^. The values chosen for the 
kinematic variables are in the range accessable at HERA. 

In order to isolate the quantitative effect of the next-to-leading correction 
to the hard scattering, in Fig. 2 we plot the same result as in Fig. 1 for the 
lowest bin along with the result obtained by including only the leading- 
order contribution to the hard scattering. The diffractive parton distributions 
and the value of the running coupling are the same in the two curves of Fig. 2. 
The effect is of the order of a factor of 2. 

Note that from the point of view of the s-channel picture of the process 
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Xp 

Figure 3: The diffractive one-jet cross section at NLO as a function of the 
hadron's momentum loss Xp for different jet transverse energies 
300 GeV; 4 GeV^ < < 70 GcV^; 0.1 < y < 0.7. Jets are defined by 
/c^-clustering, as specified in the text. 

(see, e.g., Ref. [4] and discussion in Sec. 3), the LO curve corresponds to the 
sum of the contributions 7* — > qq and 7* — > qqg, with the latter being by 
far the dominant one, since /^^^ » f^^\ The NLO curve corresponds to the 
first radiative correction to these two contributions. Fig. 2 indicates that the 
radiative corrections (besides evolution) to photon dissociation into qqg are 
important in the HERA region. 

The size of the correction raises the question of how to improve the relia- 
bility of the perturbation series. Because the short-distance matrix elements 
are weighted by parton distributions that behave very differently than the 
ordinary distributions, different physical effects at short distances will dom- 
inate the correction compared to the inclusive case. Recall from Sec. 3 that 
/^^^ stays large up to large values of (3. Then the physical cross section is 
likely to be sensitive to the behavior of the partonic cross section for small 
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Figure 4: Predictions for the jet cross section corresponding to different 
diffractivc parton distributions. The distributions in the semihard dominance 
picture (sohd curve) arc those of [8]; the HI fit distributions (dashed and 
dotted curves) are from Ref . [9] . The range of integration in and y is as 
in Fig. 3; Xp — 0.05. Inset: results for the diffractive structure function F2 
from the same sets of diffractive parton distributions. 

longitudinal momentum fraction and to the tail at finite parton transverse 

momenta. We leave the analysis of this issue to future investigation. 

In Fig. 3 we show the Xp dependence of the cross section for different 
bins. We integrate the triple-differential cross section defined above over Q^, 
with 4 GeV^ < < 70 GeV^, and plot da/ldErdxp]. Although the distri- 
butions f^^^ increase with decreasing Eq. (8), with the measured 
value a ^ 1.15 [9, 10] — , the jet cross section decreases as a result of the re- 
duction in the longitudinal phase space — see Eqs. (3) and (4). This behavior 
is due to the non-pointlike coupling of jets to the electromagnetic current, 
and characterizes the measurement of jets with respect to the measurement 
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of the structure function F2. 

Having examined the NLO perturbative corrections, let us now look at the 
impact of long-distance effects on the jet cross section. To illustrate this, we 
make use of two fits performed by the HI Collaboration to structure function 
data [9] . We compare predictions based on the diffractive parton distributions 
in the semihard dominance picture with predictions based on the diffractive 
parton distributions extracted from the HI fits. The fit distributions, while 
compatible with diffractive F2 data, are representative of different physical 
pictures of the long-distance process. 

Results are given in Fig. 4. We plot the doubly-differential jet cross 
section defined above as a function of the jet transverse energy Et- The 
solid lines in Fig. 4 correspond to the same distributions [8] as in the pre- 
vious figures; the dashed and dotted lines correspond to the fit-1 and fit-2 
distributions of HI [9]. (We recall from [9] that fit-1 distributions are quark- 
dominated, while fit-2 distributions are gluon-dominated.) The inset in the 
upper right corner of Fig. 4 shows the corresponding NLO results for the 
diffractive structure function F2, obtained from the same three sets of dis- 
tributions. The comparison gives a quantitative illustration of how much 
more sensitive the jet cross section is to long-distance effects than F2. We 
may also remark, from the results of Fig. 2 and Fig. 4, that while perturba- 
tive corrections to diffractive jet rates at HERA are large, on the order of a 
factor of 2, the effects from different scenarios for the parton distributions 
(all compatible with F2 data) are even larger, on the order of a factor of 10. 
We should likely learn a great deal about the long-distance physics of hard 
diffraction from the study of jet final states. 



5. Summary 

In this paper we have presented the factorization formula that relates 
the diffractive jet-production cross section in DIS to the diffractive parton 
distributions. Previous studies of diffractive jet leptoproduction have been 
based on approaches that do not go beyond the leading logarithm approx- 
imation. The factorization formula provides a systematic framework that 
allows arbitrarily nonleading corrections to be included. We have evaluated 
this formula explicitly to the next-to-leading order. 
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Using factorization, we have discussed how to adapt standard NLO event 
generators in order to perform NLO calculations for diffractive jet physics. 
We have used this method to compute one-jet cross sections. The method is 
general and can be applied to a variety of final-state observables. 

This improved calculational framework can be used to study the diffrac- 
tive gluon distribution. Building on previous work, and motivated by indica- 
tions from diffractive F2 data, we have discussed scenarios for long-distance 
physics in which the diffractive gluon distribution is dominated by color- 
transparency lengths of the order of the inverse of a semihard scale, Msh ~ 1 
GeV. In this case the /3 dependence becomes calculable by perturbation meth- 
ods, and leads to testable predictions for the jet cross sections. 

We recall that the factorization formula is valid up to corrections sup- 
pressed by powers of the hard scattering scale. These corrections correspond 
to multi-parton exchanges and contributions nonlinear in the parton dis- 
tributions. Since the diffractive gluon distribution is very large, this may 
overcome the power suppression. We underline the importance of searching 
for deviations from the leading power particularly in diffractive final states. 

We have observed that NLO contributions to the diffractive cross sec- 
tions are generally large in the HERA kinematic region. In the s-channel 
language in which color-dipole models are most naturally formulated, this in- 
dicates that nonleading-log corrections to photon dissociation into qqg states 
are important. We also noted the possibility that effects from finite parton 
transverse momenta in the short-distance cross section may be exposed by 
the large-/? behavior of the diffractive gluon distribution. The investigation 
of these questions is left to future work. 
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